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ABSTRACT 

(N ; 

Using observations made with the XMM-Newton Observatory, we report the 
first X-ray detection of the high magnetic field radio pulsar PSR B1916+14. We 
show that the X-ray spectrum of the pulsar can be well fitted with an absorbed 
blackbody with temperature in the range of 0.08-0.23 keV, or a neutron star 
hydrogen atmosphere model with best-fit effective temperature of ~0.10 keV, 
higher than expected from fast cooling models. The origin of the likely thermal 
emission is not well constrained by our short observation and is consistent with 
q ! initial cooling or return-current heating. We found no pulsations in these data 

£3 ! and set a la upper limit on the pulsed fraction in the 0.1-2 keV band of ~0.7. 

Implications of these results for our understanding of the different observational 
properties of isolated neutron stars are discussed. 
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ON . 1. Introduction 

O ■ 

. ^ I In the past few decades, the boom in X-ray /gamma-ray astronomy has led to the discov- 

eries of at least two new classes of isolated neutron star: magnetars and dim thermal isolated 
neutron stars (DTINS). These new classes exhibit distinctive properties different from those 
of conventional rotation-powered pulsars, the only kind of isolated neutron star known be- 
fore. Magnetars, including anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters 
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(SGRs), are slowly rotating neutron stars. They often emit luminously in the X-ray band 
and sometimes show dramatic outbursts and variability. Their X-ray emission and burst- 
ing activity are believed to be powered by enormous internal magnetic f ields (10 14 -10 15 G; 
for reviews, see IWoods fc Thompson] bood iKaspil 120071 : iMereghettil I2OO8I ) . The DTINS are 
nearby X-ray pulsars, having quasi-thermal spectra and no ra dio emission. They also appear 
to be slow rotators and are highly magnetized (B ~ 10 13 G; Kaplan fc van Kerkwiik 



2005 



Zane et al.ll2005l ; Ivan Kerkwijk fc Kaplanll2008l ; Kaplan fc van Kerkwijkll2009l ). The nature 
and origin of these DTINS are still mysteries. A third possible new class is the so-called 
"Central Compact Object" (CCOs. IPeLuca et al.l [2OO8I : IPavlov fc Lunall2009h . apparently 
isolated neutron stars located in supernova remnants. However, it is not yet clear whether 
these are a truly distinct class as some CCOs have been recently revealed to be conventional 
rotation-powered pulsars (albeit with surprisingly low B fields; e.g. the CCO 1E1207.4— 5209, 
Gotthelf fc Halpernll2007l l 



What drives the different behaviors of these neutron stars? For magnetars, the answer is 
almost certainly their enormous magnetic field. Interestingly, the DTINS also appear to have 
higher magnetic fields than those of ordinary pulsars. Hence, it is sensible to suspect that 
magnetic field could be the primary defining characteristic of the three categories of isolated 
neutron stars. Based on the study of a large sample of isolated n eutron stars , which includes 
most magnetars, some DTINS, and many ordinary radio pulsars, iPons et al.l (120071 ) reported 
an intriguing correlation between the pulsar's blackbody temperature T, determined from 
the X-ray spectrum, and the magnetic field B, inferred from the spin- down (T oc B 1 / 2 ). They 
also suggested that this correlation could be explained if the crusts of neutron stars were 
heated by magnetic-field decay, since it would significantly delay the cooling, particularly if 
the field were stronger than 10 13 G. A similar correlation was also expected if the core of the 
neutron star is heated by magnetic field decay, with heat tr a nsfer out to the surface , leading to 
an increase of the surfac e temp erature (lArras et al.l 120041 ). lAguilera et al.l (120081 ) expanded 
the work of IPons et al.l (120071 ) through 2D cooling simulations that included anisotropic 
thermal conductivity and all relevant neutrino emission processes for realistic neutron stars, 
in an attempt to unite magnetar and DTINS in a simple picture of heating by magnetic 
field. Their study shows that the DTINS could be explained if they are old neutron stars 
(~ 10 6 yr) born with 10 14 -10 15 G magnetic fields, or if they are middle-aged neutron stars 
born with 10 13 -10 14 G magnetic fields. This theory predicts that pulsars with magnetic 
field higher than 10 13 G should be hotter than is predicted by a simple cooling model with 
lower magnetic field, regardless of whether they are radio-quiet or not. Therefore, observing 
high magnetic field radio pulsars at X-ray energies and measuring the temperature of their 
thermal radiation may help us unify the different classes of isolated neutron stars. 



PSR B1916+14 is a radio pulsar having period P =1.181 s, with spin-down-inferred 
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magnetic field B = 3.2x 10 19 (PP) 1 / 2 G= 1.6 x 10 13 G, spin-down age r = P/(2P) = 8.8 x 
10 4 yr, and E = 4ir 2 IP/P' 3 = 5 x 10 33 erg s" 1 Jliulse fc Tavlor!! 19741 : Manchester et al.ll2005h . 
It is a relatively young pulsar. Given its age, PSR B1916+14 should still be hot enough to be 
X-ray de tectable, accordin g to a minimal pulsar cooling model, without magnetic-field-decay 
heating (IPage et al.ll2006l ). It is also one of the highest-magnetic- field radio pulsars known 
and may therefore be hotter because of magnetic-field decay. This makes PSR B1916+14 a 
good test subject for neutron star cooling models, and hence X-ray observations. 



2. Observations and Results 



PSR B1916+14 was obs erved by the Newton X-ray Multi-Mirror Mission (XMM ) ob- 
servatory (IJansen et al.ll200ll) on 2 008 March 25. Both the European Ph oton Imaging Cam- 
era (EPIC) pn (jStriider et al.ll200ll ) camera and the EPIC MOS cameras ([Turner et al.ll200ll ) 
were operating in full window mode with the thin filter, and with a pointing offset of 1M07. 
We analyzed the data taken in this XMM observation, and found that PSR B1916+14 was 
clearly detected in both the pn and MOS data. 



2.1. Imaging and Source Detection 

The XMM data were analyzed with the XMM Science Analysis System (SAS) version 
8.O.C0 and the latest calibrations (updated 2008 Oct 3). To exclude strong background flares 
that sometimes contaminate XMM data, we extracted light curves of photons above 10 keV 
from the entire field-of-view of the pn and MOS images, and excluded the time intervals in 
which background flares occurred for all subsequent analyses. The total exposure time of 
the observation is ~ 25 ks. However, after excluding the bad time intervals within which the 
background flux was very high (>10 counts per second) and showing significant burst-like 
features, only 12 ks of pn, 11 ks of MOS1 and 13 ks of MOS2 data were used in our analysis. 
The data were also corrected to the barycenter using the SAS barycen tool after background 
flares were exclud ed, using the nom inal pulsar position (J2000) RA 19:18:23.638(7) DEC 



+ 14:45:06.00(15) flHobbs et all 120041 ) . 



In order to find the X-ray counterpart of PSR B1916+14 , we used the SAS tool 
edetect_chain to perform a blind search for point sources. edetect_chain is designed 
to find point sources using a sliding cell technique and to calculate the significance of any 



http://xmm.esac.esa.int/sas/8.0.0/ 
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detection using a maximum likelihood method. It generates an output source list file con- 
taining information like total counts, position and significance of detected sources. In the 
pn image, a point source was detected coincident with the position of PSR B1916+14 (Fig. 
[TJ by edetect_chain. The source has 133 ± 1^| counts in 0.1-10 keV band and a likelihood 
ratio of L2 = — ln(P) ~ 143 (where P is the probability for a random Poissonian fluctuation 
to have caused the observed source counts). This source was also detected in the MOS 1 
image with 22 ± 7 counts and L 2 — 11, and in the MOS 2 images with 46 ± 9 counts and 
L 2 — 48, both in the 0.1-10 keV range. Thus, this point source was clearly detected in the 
pn and the MOS data. 

Figure [1] is the pn image, smoothed with a Gaussian profile of radius a = 8". 2. The 
small circle in the center of the image marks the radio p osition of PSR B1 916+14 . The 



position uncertainty is smaller than the size of the circle (IHobbs et al.ll2004l ). The best-fit 
position of the detected source given by edetect_chain is (J2000) RA 19:18:23.74(5) DEC 
+14:45:06.2(8), consistent with the radio position of PSR B1916+14. Therefore, it is very 
likely that the source we detected is the X-ray counterpart of PSR B1916+14 . 

The radially averaged profile of XMM 's point spread function can be approximated by 
an analytic function - the King function p(r) = A[l + (^) 2 ] _a , where p(r) is the number 
density of counts at radius r, A is a normalization parameter, r is the radial distance between 
the events and the center of the source, tq and a are parameters reflecting the size and shape 
of the point spread function (PSF) and are functions of energy and off-axis ang leE In order 
to search for evidence of extended emission from PSR B1916+14, we extracted 0.2-12 keV 
photon events from a circular region of 35" radius around the best-fit position, and calculated 
the radial distance between every photon event and the pulsar, to get the radially averaged 
profile. Using the Kolmogorov-Smirnov (K-S) test, we then compared the observed radially 
averaged profile to a model composed of a King function. Given the small off-axis angle, 
and the energy distribution of the source, we chose a = 1.6 and r = 5.25 pixels= 21". 525 
and a uniform background (0.05 photons/acrsec 2 , inferred from the 397 photons found in 
a circular background region of 50" radius and ~3' away from PSR B1916+14). The K- 
S test shows that the radially averaged profile can be well matched by the specified King 
function. Therefore, there is no evidence for extended emission near PSR B1916+14 from 
this observation. 



2 Unless otherwise specified, the uncertainties quoted in this paper represent a la range. 
3 See http://xmm.esac.esa. int/docs/documents/CAL-TN-0018- 2-6.pdf} page 6 
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2.2. Spectroscopy 

We extracted the X-ray spectrum of PSR B1916+14 from the pn data using a circular 
region of 32". 5 radius encircling the source. The source region should contain more than 
80% of the counts from a point source. The background spectrum was extracted from a 
circular region of 50" radius and ~3' away from the pulsar where no source was detected. 
Both single- and double- events were selected, but events that hit or were close to a bad pixel 
or CCD gap were excluded using the filter expression FLAG = 0&&PATTERN <= 4. A 
response file and an auxiliary response file were generated using the SAS command rmf gen 
and arf gen. The spectrum was grouped to have a minimum of 15 photons per bin using the 
f tool grppha, and was then fed to XSPEC 12.3.C0 for spectral fitting. 

We also extracted spectra from the data of the two MOS detectors using source circular 
regions of 36" radius and background regions of ~ 60" radius. Single- to quadruple- photon 
events were selected except those that landed on a bad pixel or CCD gap, using the filter 
expression of XMMEA_EM&&PATTERN <= 12. We then combined the two MOS spectra 
into a single MOS spectrum and averaged their background, response and auxiliary files using 
the f tool addspec. The resulting MOS spectrum was also grouped to have a minimum of 
15 photons per bin and was fitted jointly with the pn spectrum. 

The X-ray spectra of PSR B1916+14 can be well fit with an absorbed blackbody 
model. However, due to the small number of counts, the column density Nh was poorly 
constrained. The best-fit Njj is ~ 1 x 10 19 cm -2 , too small given the estimated distance 
and location of the pulsar. Therefore, we estimated the Njj for this pulsar based on the 
total Nh (1.58 xlO 22 cm" 2 ) of the Galaxy along the line-of-sight@ and the distance to the 
puls ar (2.1 ± 0.3 kpc, es timated from the 27.2 pc cm -3 dispersion measure of the pul- 



sar; 



Cordes fc Laziol |2002| ). and find a moderate value of ~ 0.14 x 10 22 cm -2 . Fixing Nh 
to this value, the best-fit blackbody temperature for the 0.1-2 keV spectra (Fig. [2]) is 
0.13 ± 0.01 keV (Table [p. The model-predicted absorbed flux in the 0.1-2 keV range is 
1.4±0.3xl0~ 14 ergs s _1 cm -2 . Assuming a distance of 2.1 kpc, we find the bolometric X-ray 
luminosity of PSR B1916+14 to be ~ 3 x 10 31 ergs s -1 . 

By fixing Nh while fitting the pn and MOS spectra, we underestimate the uncertainties 
of the best-fit parameters. In order to get a sense of the real uncertainty of kT, we tried fitting 
with a range of N H . N H likely lies between 0.07x 10 22 cm -2 and 0.3 x 10 22 cm -2 . It is probably 
not smaller than 0.07 x 10 22 cm -2 because the distance estimated from the dispersion measure 



http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
5 http:// cxc.harvard.edu/toolkit / colden.jsp 
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is unlikely to be incorrect by more than 50%. Also, an absorbed blackbody model with Nh 
higher than 0.3 x 10 22 cm -2 cannot fit the spectra well for any kT. With Nh restricted to 
lie between these two values, the acceptable (null hypothesis possibility of the fit > 0.01) 
range of kT is 0.08 to 0.23 keV (with blackbody radius range from ~6 km to ~0.2 km). This 
temperature range reflects reasonable uncertainties on Nh, and is quoted in the abstract and 
Figures [3] and 0] (see below) . 

The pn and MOS spectra could also be fit with a power-law model, with a best-fit Nh 
of 0.12+°;^ x 10 22 cm" 2 and a photon index of T ~3.5±J;f (Tabled]). The lack of source 
photons with energy above 2 keV results in a soft best-fit power-law model. This is rarely 
seen from other non-thermally emitting pulsars. Therefore, it is very unlikely that the X-ray 
emission of PSR B1916+14 is non-thermal. 

A neu tron-star hydrogen atmosphere (NSA) model (with magnetic field strength set 



to 10 13 G; IZavlin et al.l Il996l ; iPavlov et al.l Il995l ) could also fit the pn and MOS spectra. 



However, the parameters are even less constrained in comparison with the blackbody model. 
We had to freeze the mass of the neutron star to 1.4 M & and the distance to 2.1 kpc to get 
a better-constrained fit (xV) = 14.3(17)). The best-fit N H is 0.23ig;{g x 10 22 cm" 2 , kT is 
0.10 ± 0.04 keV, and the resulting best-fit neutron-star radius is ~6 km. Unfortunately, the 
radius is not well constrained in this model and has a la upper limit of 20 km. The model 
predicted 0.1-2 keV X-ray luminosity is (7 ± 4) x 10 31 ergs s -1 . 



2.3. Timing analysis 

To search for X-ray pulsations from PSR B1916+14 , we folded all the pn source events 
from a total 25 ks exposure without filtering for the background flares using 16 phase bins and 
a contemporaneous ephemeris which was derived from radio timing data obtained using the 



76-m telescope at the Jodrell Bank Observatory (IHobbs et al.ll2004l ). The MOS full- window 
mode data is useless for timing analysis because of its 2.7 s time resolution. A total of 374 
pn photon events were extracted without filtering for background flares, all in the range of 
0.1-2 keV, from a source region of 15" radius chosen to reduce the number of background 
photons and improve the signal-to-noise ratio. From the same energy band and the same 
CCD, 1945 events were found in a circular background region of 50" radius where no source 
was detected by the SAS tool edetect_chain. If the background is uniformly distributed, 
there should be 175±4 background photons in the source region. The folded and binned 
light curve was fit to a constant line. The best-fit \ 2 was 9-6 for 15 degrees of freedom. 
Therefore no significant pulsations were detected. 
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In order to determine an upper limit on the pulsed fraction, we simulated event lists 
with the same total number of counts as in the observed event list. The simulated event 
lists were generated assuming the signal has a sinusoidal profile starting at a random phase 
and a specified area pulsed fraction, where the area pulsed fraction is defined as the ratio 
of the pulsed part of profile to the entire profile. For a sinusoidal profile F = A sin 8 + B, 
where F is the count rate, A is the amplitude, B is the DC level, and 9 is the phase, the area 
pulsed fraction is A/(2B + A). By specifying different area pulsed fractions, we found that, 
if we set the area pulsed fraction of the simulated event lists to 35%, then approximately 
68% of them would be detected with > 3a significance. Because there are ~175 background 
photons in the 374 photons from the source region, the la area pulsed fraction upper limit 
of the pulsar is ~ 0.7. This is not an especially interesting constraint because the number 
of source photons was so small that even a highly pulsed signal could have gone undetected. 



Discussion 



The X-ray spectrum of PSR B1916+14 is soft, and therefore is most likely thermal. 
The blackbody temperature lies in the range of 0.08-0.23 keV, and the best-fit effective 
temperature of NSA model is 0.10±0.04 keV. Given PSR B1916+14 's age, its temperature is 
consistent with what one would expect (kT ~ Q. 07-0. 11 keV) from minimal cooling models in 
which magnetic field is not considered (IPage et al.ll2006l ). Its estimated bolometric luminosity 
(~ 3 x 10 31 erg s _1 ; Table [1]) is somewhat low when compared with the curves of Page et 
al. (2006), suggesting that the pulsar may have a light-element envelope. However the 
substantial uncertainties on the luminosity preclude a firm conclusion. On the other hand, 
fast cooling models predict much lower temperatures (kT < 30 eV; e.g. Yakovlev & Pethick 
2004) which would be undetectable with current instruments. Therefore, fast cooling seems 
unlikely for PSR B1916+14 from our observation. 

The small best-fit blackbody radius (0.8 ±0.1 km; Tabled]) suggests polar-cap reheating. 
If the emission of the pulsar is due to c urvature radiation, return-c urrent heating is predicted 
to giv e rise to an X-ray luminosity (see lHarding fc Muslimovll200ll and Eq. 7.2 in lKaspi et al. 
20061 for details) 



L\ ' ~ 10 erg s 



0.4p-6/7( T /l 6)-lA 

l.OP-Va 



ifP < 0.1(P/10 12 ) 4 / 9 
ifP > 0.1(P/10 12 ) 4 / 9 , 



where r is in yr, P is in s, and B is in G. For PSR B1916+14 , l}+ R ^ is ~ 9 x 10 30 erg s _1 , only a 
factor of ~3 smaller than the estimated blackbody bolometric luminosity (see Table [1]). If the 
emission is due to inverse Compton scattering, the return current heating will be much less 
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effective ( IHarding &: Muslimovll2002l ). Given that the best-fit blackbody parameters are not 
well constrained, we cannot rule out return current heating as the origin of PSR B1916+14 's 
X-ray luminosity. We note that the NSA model (Section 12.21) yields a larger radius, although 
it is also not well constrained. 

To compare the properties of PSR B1916+14 to those of other X-ray-detected radio 
pulsars, we have collected the temperature, magnetic field strength, and spin-down energy 
of a dozen such pulsars from the literature, and listed them in Table [2J For pulsars in this 
Table, we also made plots of their temperature versus age and magnetic field (Figs. [3] and HJ 
respectively). Given the large uncertainty on the temperature measurement from our short- 
exposure observation of PSR B1916+14 , we cannot conclude here whether thermal emission 
is consistent with minimal cooling or if the neutron star is hotter than lower-magnetic-fiel d 
pulsars of the same age, as expected in some models (IPons et al.ll2007l ; lAguilera et al.ll2008l ). 
A longer observation in the future may be able to distinguish among thermal models. 

However, we note with interest that the previously published temperature of PSR 
J0538+2817 is surprisingly high, in spite of its relatively large age (40 kyr) and relatively low 
magnetic field (7.3 x 10 11 G; Table EJ. In contrast to PSR B1916+14, PSR J0538+2817's 
emission is unlikely to be from polar-cap reheating because of the very high required ef- 
ficiency for conversion of E to X-ray luminosity, ~ 10~ 2 , compared with the ~ 5 x 10~ 4 
predicted by Equation 1 above for this pulsar. If correct, the high temperature suggests a 
wider range of possible tempera tures for young; ne utron stars than is currently predicted. 
This would be a challenge to the IPons et al.l (120071 ) model. 



Although PSR B1916+14 is detected by XMM as a point source with no evidence of 
extended emission, it is still possible that extended emission was too faint to be detectable. 
Based upon the number of counts in the background region, we found that, in the 0.3-8 keV 
band, extended emission of surface brightness smaller than ~ 3 x 10~ 6 count s _1 arcsec -2 
would not be detected (with 3a significance) in our observation. This limits our sensitivity 
for detecting a very faint pulsar wind nebula (PWN) like th at of the Geminga pulsar, which 



has a surface brightness of ~ 1 x 10 6 count s 1 arcsec 2 (jPavlov et al.l 120061 ) in the same 



energy range, despite the fact that Geminga is closer by a factor of ~8. Assuming there is 
an undetected PWN around PSR B1916+14 having a spectrum like that of the Geminga 
PWN (power law with index 1.0), we can estimate the upper limit of its surface brightness 
to be ~ 3 x 10~ 17 erg cm _2 s _1 arcsec -2 , in the 0.3-8 keV range. If we further assume that the 
PWN is uniformly distributed in a circular region of radius 20", then this surface brightness 
upper limit corresponds to a PWN luminosity upper limit of ~3xl0 31 erg s . 



Geppert et al.l (120041 ) showed that the presence of a very high magnetic field could cause 



inhomogeneous thermal conductivity in the neutron star crust and lead to the formation of 
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hot spots on the neutron star surface. It could also cause strong radiative beaming of the 
thermal emission from the neutron star. These effects could give rise to highly pulsed X- 
rays, as in the 74 ± 14% pulsed fraction observed from the high-magnetic- field radio pulsar 



PSR J1119— 6127 ([Gonzalez et al.l 120051 ). However, limited by the exposure time of the 
observation, we did not detect any X-ray pulsations from PSR B1916+14 . Future longer 
observations will be useful for better constraining its pulsed fraction. 

We thank S. Bogdanov and S. Guillot for useful comments. We thank the referee Craig 
Heinke for helpful suggestions. VMK is supported by an NSERC Discovery Grant, CIFAR, 
FQRNT, the Canada Research Chairs Program, and the McGill Lome Trottier Chair. 
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Table 1. Best-fit Spectral Parameters for PSR B1916+14 



Parameter Blackbody modeP Power-law model NSA model 



N H (10 22 cm~ 


2 ) 


0.14 


12+ 005 


n O o+0.09 
u - ZJ -0.04 


kT (keV) 




0.13±0.01 




0.10 ±0.04 


Rbb (km) 




0.8 ±0.1 






R NS (km) 








~6 


r 






o r+1.6 
°-°-0.7 




X 2 (dof) 




14.1(18) 


13.9(17) 


14.3(17) 


/abs b (ergs s _1 


cm^ 2 ) 


1.4±0.3xl0- 14 


1.7^; 4 xl0- 14 


(1.7±0.3) x 10- 14d 


/unabs c (ergs s 


- 1 cm -2 ) 


(5±l)xl0- 14 


- 2xl0- 13 


(1.4±0.8) x 10- 13d 


L x (ergs s _1 ) 




- 3xl0 31c 


- Ixl0 32f 


(7 ±4) x 10 31d 



a Best-fit parameters of absorbed blackbody fit to the XMM spectra. Njj was frozen 
when fitting, so the uncertainties of the parameters, especially that of the kT, do not 
reflect the uncertainties on Nh] see text for details. Emission radius Rbb was inferred 



assum ing a distance of 2.1 kpc (estimated from the dispersion measure; iHobbs et al 



2004). 



b Absorbed X-ray flux, / a b s in the 0.1-2 keV range. 

c Unabsorbed X-ray flux, / un abs, in the 0.1-2 keV range; the uncertainty was propagated 
from the uncertainties on the parameters and absorbed flux. 

d When fitting with NSA model, the flux and luminosity are estimated using the cf lux 
model in xspec. 

°Bolomctric X-ray luminosity derived assuming a distance of 2.1 kpc. 

f Inferred X-ray luminosity in the 0.1-10 keV range assuming a distance of 2.1 kpc. 



Table 2. Parameters of the X-ray-detected Radio Pulsars 



Name 


Age a (kyr) 


kT h (eV) 


B c (G) 


E c (ergs 4 ) 


Observatory 


Spectral model d 


Reference 


B1916+14 


88 


135(45) 


1.6xl0 13 


5.1 x 10 33 


XMM 


BB 


this work 


B1055-52 


540 


68(3) 


l.lxlO 12 


3.0 x 10 34 


XMM 


BB+BB+PL 


De Luca et al. 2005 


J0633+1746 (Geminga) 


340 


41.4(0.1) 


1.6xl0 12 


3.2 x 10 34 


XMM 


BB+PL 


Jackson & Halpcrn 2005 


B0656+14 


110 


56.0(0.9) 


4.7xl0 12 


3.8 x 10 J 


XMM 


BB+BB+PL 


De Luca et al. 2005 


B0355+54 


562 


82(4) 


8.4xl0 n 


4.5 x 10 34 


RO SAT /Einstein 


BB 


Slane 1994 


J0538+2817 


40 c 


183(3) 


7.3x10" 


4.9 x 10 34 


XMM 


BB 


McGowan et al. 2003 


B2334+61 


41 


109(35) 


9.9xl0 12 


6.2 x 10 34 


XMM 


BB+PL 


McGowan et al. 2006 


B1823-13 


21 


135(14) 


2.3xl0 12 


2.8 x 10 36 


Chandra 


BB 


Pavlov et al. 2008 


B1706-44 


17.4 


143(14) 


3.1xl0 12 


3.4 x 10 36 


Chandra 


BB+PL 


Gotthelf et al. 2002 


J1811-1925 


2 f 


<150 


1.7xl0 12 


6.4 x 10 36 


Chandra 


BB+PL 


KasDi et al. 2006 


B0833-45 (Vela) 


11 


91(3) 


3.4xl0 12 


6.9 x 10 36 


XMM 


BB+BB+PL 


Manzali et al. 2007 


J0205+6449 


2.4 f 


112(9) 


3.6xl0 12 


2.7 x lO 37 


Chandra 


BB+PL 


Slane et al. 2004 
WcisskoDf et al. 2004 


B0531+21 (Crab) 


0.955 f 


<172 


3.8xl0 12 


4.6 x lO 38 


Chandra 


BB+PL 



a The spin-down age unless otherwise noted. 

b The blackbody temperature or the temperature of the softer blackbody component as measured by fitting the data with different spectral 
models as listed in this Table. 



c Numbers were found in the ATNF database ((Manchester et al.ll2005l ). 

d BB: blackbody model; BB+PL: blackbody plus power-law model; BB+BB+PL: two blackbody plus power-law model. 



3 The age of PSR J0538+2817 estimated based on the proper motion of the pulsar from its associated supernova remnant (SNR) (|Ng et al 



2007J)- 



f The age of SNR with wh ich the pulsar is associated, estimated based on its expansion rate: PSR J1811— 1925 ( Kaspi et al.l I2006T ) : PSR 
J0205+6449 |Chevalierll20oi) . The age of PSR B0531+21 (Crab pulsar) is based on historical record. 
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Fig. 1. — XMM image in the 0.2-2 keV energy band, smoothed by a Gaussian profile of 
8.2" radius (a profile that is slightly oversampling the telescope's PSF). The radio position 
of PSR B1916+14 is labeled by the white circle. Note that the radius of the circle is much 
larger than the uncertainty on the radio position. 
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Fig. 2. — XMM spectra (upper is pn, lower is combined MOS) of PSR B1916+14, with the 
best blackbody fit (see Tabled]). The spectra are binned to contain a minimum of 15 counts 
per bin. 
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Fig. 3. — Observed temperature (kT) versus age for X-ray-detected radio pulsars. Note 
that for PSR B 19 16+ 14 we used the kT measured by allowing iV# to vary in a reasonable 
range; the same is not necessarily true for the other measurements; see original references 
for details. 
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Fig. 4.— Observed temperature (kT) versus magnetic field strength for X-ray-detected 
radio pulsars. See caption for Figure [3] for caveats. 



